THE KISSING PROBLEM IN THREE 
DIMENSIONS 



Oleg R. Musin * 



Abstract 

The kissing number k (3) is the maximal number of equal size nonover- 
lapping spheres in three dimensions that can touch another sphere of the 
same size. This number was the subject of a famous discussion between 
Isaac Newton and David Gregory in 1694. The first proof that fc(3) = 12 
was given by Schiitte and van der Waerden only in 1953. In this paper 
we present a new solution of the Newton-Gregory problem that uses our 
extension of the Delsarte method. This proof relies on basic calculus and 
simple spherical geometry. 

Keywords: Kissing numbers, thirteen spheres problem, Newton-Gregory 
problem, Legendre polynomials, Delsarte's method 

1 Introduction 

The kissing number k(d) is the highest number of equal nonoverlapping spheres 
in R d that can touch another sphere of the same size. In three dimensions the 
kissing number problem is asking how many white billiard balls can kiss (touch) 
a black ball. 

The most symmetrical configuration, 12 billiard balls around another, is 
achieved if the 12 balls are placed at positions corresponding to the vertices of 
a regular icosahedron concentric with the central ball. However, these 12 outer 
balls do not kiss each other and may all be moved freely. So perhaps if you 
moved all of them to one side a 13th ball would possibly fit in? 

This problem was the subject of a famous discussion between Isaac Newton 
and David Gregory in 1694 (May 4, 1694; see interesting article ;21 for details 
of this discussion). Most reports say that Newton believed the answer was 12 
balls, while Gregory thought that 13 might be possible. However, Casselman 
jS] found some puzzling features in this story. 

This problem is often called the thirteen spheres problem. Hoppe [5] thought 
he had solved the problem in 1874. But there was a mistake - an analysis of this 
mistake was published by Hales in 1994 |S] (see also [201 ) - Finally this problem 
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was solved by Schiitte and van der Waerden in 1953 ^Hj. A subsequent two- 
page sketch of an elegant proof was given by Leech JJJ in 1956. Most people 
agree that Leech's proof is correct, but there are gaps in his exposition, many 
involving sophisticated spherical trigonometry. (Leech's proof was presented in 
the first edition of the well-known book by Aigner and Ziegler pQ, the authors 
removed this chapter from the second edition because a complete proof would 
have had to include so much spherical trigonometry.) The thirteen spheres 
problem continues to be of interest, and new proofs have been published in the 
last few years by Hsiang ^U], Maehara Boroczky [3] and Anstreicher PJ. 

The main progress in the kissing number problem in high dimensions was at 
the end of the 1970s. Levenshtein an d independently, Odlyzko and Sloane 
[TBI (= Chap.13]) using Delsarte's method in 1979 proved that fc(8) = 240 
and fc(24) = 196560. This proof is surprisingly short, clean, and technically 
easier than all proofs in three dimensions. However, d — 8, 24 are the only 
dimensions in which this method gives a precise result. For other cases (for 
instance, d = 3,4) the upper bounds exceed the lower. 

We found an extension of the Delsarte method in 2003 ^l](see details in 
that allowed us to prove the bound fc(4) < 25, i.e. fc(4) = 24. This extension 
also yields a proof k(3) < 13. 

The first version of these proofs was relatively short, but used a numerical 
solution of some nonconvex constrained optimization problems. Later ^5] these 
calculations were reduced to calculations of roots of polynomials in one variable. 

In this paper we present a new proof of the Newton-Gregory problem. This 
proof needs just basic calculus and simple spherical geometry. 

2 fc(3) = 12 

Let us recall the definition of Legendre polynomials Pk{t) by the recurrence 
formula: 

3,1 2k- 1 k- 1 
Pa = 1, P 1= t, P 2 = -t 2 - P k = —j—tPk-! — P k -2 

or equivalently 

Pk(t) = ^ ^(t 2 - l) fe (Rodrigues' formula). 

Lemma 1. Let X = {x%, x%, . . . , x n } be any finite subset of the unit sphere S 2 
in R 3 . By <pij = dist(cci , Xj) we denote the spherical (angular) distance between 
Xi and Xj . Then 

n n 



2 



This lemma easily follows from Schoenberg's theorem f° r Gegenbauer 
(ultraspherical) polynomials • (Note that P^ = G^' .) For completeness we 
give a proof of Lemma 1 in the Appendix. 

Let 

,. , 2431 o 1287 7 18333 - 343 . 83 , 213 2 t 1 

fit) = t 9 t 7 + 1 5 + 1 4 t 3 t 2 + . 

J w 80 20 400 40 10 100 10 200 

Remark. This polynomial of degree 9 is satisfying the assumptions of the ex- 
tended Delsarte's method An algorithm for calculating suitable poly- 
nomials is presented in the Appendix of |15|. 

Lemma 2. Suppose X — {xi,X2, ■ ■ ■ ,x n } C S 2 . Then 

n n 

s ( x ) -EE^ cos (<M) ^™ 2 - 

i=i j=i 

Proof. The expansion of / in terms of Pk is 

/ - ± c t P„ . P„ + |P, + |« + |fl, + §P 4 + + |p 8 . 
We have Co = 1, Cfc > 0, fc = 1, 2, . . . , 9. Using Lemma 1 we get 

9 n n n n 

s{x) = £ c *EE iMco S (^)) ^ E E c ° p o = 

fe=0 i=l j'=l i=l j=l 

□ 

If n unit spheres kiss the unit sphere in R 3 , then the set of kissing points is an 
arrangement on the central sphere such that the (Euclidean) distance between 
any two points is at least 1. So the kissing number problem can be stated in 
another way: How many points can be placed on the surface of S 2 so that the 
angular separation between any two points is at least 60°? 

Lemma 3. Suppose X — {x\, X2, ■ ■ ■ , x n } is a subset of S 2 such that the angular 
separation (f>ij between any two distinct points Xi,Xj is at least 60°. Then 

n n 
i=l 3=1 

We give a proof of Lemma 3 in the next section. 
Theorem. fc(3) = 12. 

Proof. Suppose A is a kissing arrangement on S 2 with n — fc(3). Then A 
satisfies the assumptions in Lemmas 2 and 3. Therefore, n 2 ^ S(X) < l'Sn. 
From this n < 13 follows, i.e. n ^ 12. From the other side we have fc(3) ^ 12, 
showing that n = fc(3) = 12. □ 
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3 Proof of Lemma 3. 

We need one fact from spherical trigonometry, namely the law of cosines: 

cos 4> — cos Q\ cos 9 2 + sin 9\ sin 9 2 cos ip, 

for a spherical triangle ABC with sides of angular lengths 6\, 9 2 , 4> an d 
ZBAC = ip (Fig. 1). For ip = 90°, this reduces to the spherical Pythagorean 
theorem: cos <p — cos 0\ cos 9 2 ■ 

Proof. 1. The polynomial f(t) satisfies the following properties (see Fig. 2): 
(i) f(t) is a monotone decreasing function on the interval [— 1, —to}; 

/(*) < for te (-to, 1/2]; 
where f(-t ) = 0, t w 0.5907. 

These properties hold because f(t) has only one root — t on [—1, 1/2], and 
there are no zeros of the derivative f'(t) (eighth degree polynomial) on [—1, —to]- 




Let Si(X) := ^/(cos(c^)), then = ^^(X). From this it follows 

j=i i=i 



that if Si(X) < 13 for i = 1, 2, . . . , n, then < 13n. 

We obviously have 0^ = 0, so /(cos^j) = /(l). Note that our assumption 
on X ((f)ij > 60°, i 7^ j) yields cos^j < 1/2. Therefore, cosfaj lies in the 
interval [-1,1/2]. By (ii) we have /(cos^j) ^ whenever cos 4>ij € [—to, 1/2]. 
Let J(i) := {j : coscpij e [—1, —to)}- We obtain 

< T,(X) := /(l) + J] /(cos&j). (1) 

Let 6» = arccosto w 53.794°. Then j e J(i) iff 4J > 180° - 6 , i.e. 9 3 < 6 , 
where 9j = 180° — fcj. In other words all Xij, j G J(i), lie inside the spherical 
cap of center eo and radius 0q, where eo = — x, is the antipodal point to £j. 
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2. Let us consider on S 2 points eo, y\, . . . , y m such that 

(f) itj := dist(j/i, j/j) ^60°, V i ^ j, 0j := dist(e , y») < 9 for 1 < i < to. (2) 

Denote by /x the highest value of m such that the constraints in (2) allow a 
nonempty set of points yi, . . . , y m . 

Suppose that < m ^ [i and Y = {t/i, . . . , y m } satisfies (2). Let 

H{Y) = H( Vl , . . . , y m ) := /(l) + /(- cosfc) + . . . + /(- cos0 m ), 

h m := sup{iJ(Y)}, h max :=max{/i ,/ii,...,/i^}. 
y 

It is clear that Ti(X) ^ /i m , where m = |J(i)|. From (1) it follows that 
Si(X) < /i TO . Thus, if we prove that h max < 13, then we prove Lemma 3. 

3. Now we prove that fi ^ 4. 

Suppose Y = {yi, . . . ,y m } C S 2 satisfies (2). By symmetry we may assume 
that eo is the North pole and y, has polar coordinates (9i,<pi). Then from the 
law of cosines we have: 

cos fcj = cos $i cos 0j + sin 9i sin 0j cos(<y9i — <fij). 

Note that 0, > for m ^ 2. Conversely, t/j = eo, 9j — faj ^ 60° > 6o, a 
contradiction. From (2) we have coscpij < 1/2, then 

. . 1/2 — cos 6i cos ,„. 

cos(<^ -<fj)<: . . -. 3) 

sin t'i sin 6*j 

Let 

„ s 1/2 — cosacos/3 

:= : — a . 

sin a sm p 

then 

. dQ(a,[3) 2 cos f3 — cos a 

Q a ( a iP) = — o = "7T~r^ — j-- 

act 2 sin a sm /3 

From this it follows that if < a,/3 < 6*o, then cos/3 > 1/2 (because 9q < 60°); 
so then Q' a (a, f3) > 0, and Q{a, 0) < Q(0 O , /?) = Q(0, 9 ) < Q(0 O , 0o). Therefore, 

1/2 - cos 6»i cos 6j 1/2 - cos 2 O _ 1/2 - 
sin #j sin 0j sin 2 # 1 _ *o 

Combining this inequality and (3), we get 

COS{(pi -ipj)^ — 7^. 

i — r 

Note that arccos((l/2 - - to)) « 76.582° > 72°. This implies that m sC 4 

because no more than four points can lie in a circle with minimum angular 
separation between any two points greater than 72°. 
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4. Now we have to prove that h max = max{/io, hi, h 2 , ft.3, /14} < 13. 
We obviously have h = /(l) = 10.11 < 13. 

From (i) follows that f(-cosO) is a monotone decreasing function in 9 on 
[0, 80}. Then for m = 1 : H{y\) = /(l) + /(— cos#i) attains its maximum at 
6>i = 0. So then 

fci = /(!) + /(-!) = 12-88 < 13. 

5. Let us consider for m = 2, 3, 4 an arrangement {eo, j/i, . . . , j/ m } in S 2 that 
gives H(Y) = h m . Here ?/i 7^ eo (see 3). Note that in this arrangement, points 
yk cannot be shifted towards eo because in this case H{Y) increases. 

For tyl — 2 this yields eo <E J/1J/2, and dist^i, 2/2) = 60°. If e ^ 2/1Z/2, then 
the whole arc 2/12/2 can be shifted towards eo- If dist(j/i, y 2 ) > 60°, then j/i (and 
2/2) can be shifted towards eo- 

For m = 3 we prove that A3 = 2/12/22/3 is a spherical regular triangle with edge 
length 60°. As above, eo € A3, otherwise the whole triangle can be shifted to- 
wards eo- Suppose dist(j/i, yi) > 60°, i = 2, 3, then dist(t/i, eo) can be decreased. 
From this follows that for any j/j at least one of the distances {dist(j/j, yj)} is 
equal to 60°. Therefore, at least two sides of A3 (say 2/12/2 and 2/12/3) have length 
60°. Also dist(2/2, 2/3) = 60°, conversely 2/3 (or 2/2, if e € 2/12/3) can be rotated 
about 2/1 by a small angle towards e (Fig. 3). 

For Tji — 4 we first prove that A 4 := conv Y (the spherical convex hull of Y) 
is a convex quadrilateral. Conversely, we may assume that 2/4 <E 2/12/22/3- 

The great circle through 2/4 that is orthogonal to the arc e 2/4 divides S 2 into 
two hemispheres: Hi and Hi- Suppose eo € Hi, then at least one 2/j (say 2/3) 
belongs to H 2 (Fig. 4). So the angle Zeo2/42/3 greater than 90°, then (again 
from the law of cosines) dist(j/3, eo) > dist(j/3, 2/4). Thus, 
63 = dist(2/3,eo) > dist(?/3, 2/4) 5^ 60° > 60 — a contradiction. 

Arguing as for m = 3 it is easy to prove that for any vertex yi there are at 
least two vertices yj at the distance 60° from 2/i- Note that the diagonals of A4 
cannot be both of lengths 60°. Conversely, at least one side of A 4 is of length 
less than 60°. Thus, A 4 is a spherical equilateral quadrangle (rhomb) with edge 
length 60°. 

6. Now we introduce the function Filtp), 1 where ip € [60°, 29o]- 

iq(V) := max {Fi(M)}, = /(- cos0) + /(- costy - 0)). 

So if dist(y i,yj) = ip, then 

/(-cosfc) + /(-cos0 j ) ^ Fi(V). (4) 

Therefore, 

H(yi,y 2 ) ^h 2 = f(l) + Fi(60°) w 12.8749 < 13. 

1 For given i/>, the value Fi{ip) can be found as the maximum of the 9th degree polynomial 
(l(s) = Fi(6,4>), s = cos (6 - if)/2), on the interval [cos(6» - tp/2), 1]. 
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7. Let m = 4, di = dist(j/i,t/ 3 ), d 2 = dist(y 2 , 2/4), Since A 4 = 2/12/22/32/4 is 
a spherical rhomb, we have cos(di/2) cos(d 2 /2) = 1/2 (Pythagorean theorem, 
the diagonals 2/12/3, 2/22/4 01 A4 are orthogonal). So if 



p(s) := 2 arccos 



1 



2cos(s/2)' 



then 



p(d 1 ) = d 2 , P (d 2 ) = d 1 , p(90°) = 90°, /9(p(s))=s. 
Suppose d\ ^ d 2 . The inequalities #j $5 #0 yield d 2 ^ 2#o- Then 

p(26»o) < di 90° rf 2 s$ 26» - 

Now we consider two cases: 
1) p(26> ) sC di < 77°, and 2) 77° d 1 < 90°. 

1) Clearly, Fi(ip) is a monotone decreasing function in tp. Then (4) implies 

f(-coB0 1 ) + f(-coa0 3 )^F 1 (d 1 )^F 1 (p(26 o )), 
/(-cos0 2 ) + /(-cos0 4 ) < fi(d 2 ) = Fi(p(di)) < Fi(p(77°)), 

so then 

if(F) < /(l) + F 1 (p(2e )) + F!(p(77°)) « 12.9171 < 13. 

2) In this case we have 

H(Y) ^ /(l) + Fi(77°) + Fi(90°) w 12.9182 < 13. 
Thus, /14 < 13. 

8. Our last step is to show that h 3 < 13. 2 

Since A 3 is a regular triangle, H(Y) = /(l) + /(— cos#i) + /(— cos# 2 ) + 
/(—cos 6*3) is a symmetric function in the $i, so it is sufficient to consider the 
case 6*i < 62 < #3 ^5 #o- 

In this case i?o ^ $3 ^ #0, where R = arccos y/2/3 ~ 35.2644° is the 
(spherical) circumradius of A3. 

2 A more detailed analysis shows /13 12.8721, /14 » 12.4849. 
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Let y c be the center of A3. We have 7 := Zyiy 3 y c = /Ly2y^y c - Using the low 
of cosines for the triangle y\y^y c , we get 7 = arccos y2/3, i.e. j = Rq. 
Denote the angle Zeoysy c by u. Then (see Fig. 5) 

cos 61 = cos 60° cos 63 + sin 60° sin #3 cos (Rq — u), 

cos 62 = cos 60° cos 63 + sin 60° sin #3 cos (Rq + u) . 

where ^ u ^ uq := arccos(cot 63/VS) — i?o- Note that if u = u , then 
8 2 = 3 ; u = yields 9 X = #2; and if < u < u , then 9 1 < 9 2 < 9 3 . 

For fixed 63 = ip, H(y±, y 2 ) is a polynomial of degree 9 in s = cosu. Denote 
by -FMVO the maximum of this polynomial on the interval [cosuo, 1]. 

Let 

{Vi, • . . , Ve} = {Ro, 38°, 41°, 44°, 48°, 9 }. 

It is clear that F 2 (ip) is a monotone increasing function in ip on [Rq, 9q\. From 
the other side, /(— cos ip) is a monotone decreasing function in ip. Therefore for 
6*3 e [ipi,ipi+x] we have 

H(Y) = H(yi,y 2 ) + f(-cos9 3 ) < Wl := F 2 (V> i+ i) + /(- cos Vi). 

Since, 

{ioi, . . . , w 5 } w {12.9425, 12.9648, 12.9508, 12.9606, 12.9519}, 

we get /13 < max{wj < 13. 

Thus, h m < 13 for all m as required. □ 

Appendix. Proof of Lemma 1. 

In this proof we are using Schoenberg's original proof [T5j which is based on 
the addition theorem for Gegenbauer polynomials. 3 The addition theorem for 
Legendre polynomials was discovered by Laplace and Legendre in 1782-1785: 

Pfe(cos 9\ cos 9 2 + sin 9\ sin 92 cos <p) 

= P k (cosd 1 )P k (cose 2 ) + 2 V ~ m i; Pr(cosW m (cosfl 2 ) costo^ 

* — ' (fc + to)! 

m— 1 v 7 

fc 

= 51 Crn^Pk" (COS 6-i)P^ (COS 9 2 ) COSmip, 
m=0 

where 

pn*) = (i-* 2 )^^ fc (*)- 

(See details in g] and 0.) 

3 Pfender and Ziegler [T7l give a proof as a simple consequence of the addition theorem for 
spherical harmonics. This theorem is not so elementary. The addition theorem for Legendre 
polynomials can be proven by elementary algebraic calculations. 
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Proof. Let X = {x\, . . . , x n } C S 2 and Xi has spherical (polar) coordinates 
[9i, ifi). Then from the law of cosines we have: 

cos <j)ij — cos Oi cos 8j + sin 9i sin 9j cos ifiij , ifiij :— ifi — ipj , 

which yields 

k 

^P fe (cos(/) 4J ) = c m .kPk m (cos9i)Pk l ( cos0 3) cosm <fij 

i,j i.j m—0 

= ^ C mM ^ Um,iUm,j COSmipij, U mt i = P™(cOS#;). 

m i.j 

Let us prove that for any real u\, . . . , u n 

UiUj cosmipi.j ^ 0. 

Pick n vectors v\, . . . ,v n in R 2 with coordinates Vi = (cosrn^,, sinm<pi). If 
v = uivi + . . . + u n v ni then 

^ IM| 2 = (v,v) = UjUj cos rrnpij . 
This inequality and the inequalities c m .k > complete our proof. □ 
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